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Abstract 


We have analyzed a large number of hard X-ray events which were recorded by the 
Hard X-Ray Burst Spectrometer (HXRBS) on the Solar Maximum Mission (SMM) during 
the maximum of the 21 st solar cycle (circa 1980) in order to study their statistical correla- 
tion with type III bursts. We confirm qualitatively the earlier finding by Kane (1981) that 
flares with stronger hard X-ray emission, especially those with harder spectra, are more 
likely to produce a type III burst. The observed distribution of hard X-ray and type III 
events and their correlations are shown to be satisfactorily described by a bivariate distri- 
bution consistent with the assumption of statistical linear dependence of X-ray and radio 
burst intensities. From this analysis we determine the distribution of the ratio of X-ray 
intensity (in counts/s) to type III intensity (in solar flux units) which has a wide range 
and a typical value for this ratio of about 10. The implications of our results for impulsive 
phase models are discussed. 

Subject headings: Sun: flares - Sun: radio radiation - Sun: x-rays 
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I. Introduction 


It is generally accepted that the electrons which produce the impulsive hard X-ray 
bursts and those which cause the type III radio emissions during a solar flare have a 
common origin. The hard X-rays are a result of bremsstrahlung radiation of electrons 
which penetrate to the lower corona and chromosphere while type III radiation is produced 
by those electrons which escape on open field lines generating plasma turbulence as they 
travel away from the sun. The correlation of hard X-ray emission with type III emission has 
been studied both statistically (Kane 1981) and by detailed analysis of individual events 
(see e.g. Kane, Pick, and Raoult 1980, Kane, Benz, and Treumann 1982, Dennis et al. 1984, 
or Raoult et al. 1985). There is a wide variation in the relative strengths of the hard X-ray 
and type III bursts. In fact, the majority of hard X-ray events are not accompanied by an 
identifiable type III burst and vice versa. From this data one can estimate the ratio of the 
number of electrons involved in the production of hard X-rays to those producing type III 
radiation. For example, Ramaty et a/. (1980) estimates an average value of 10 2 to 10 4 for 
this ratio. However, the range of this ratio is very large. In some hard X-ray events with 
no detectable type III burst this number must exceed 10 7 (Simnett and Benz 1986). 

This kind of information is not very useful in the theoretical understanding of the 
phenomena. A knowledge of the mean value and dispersion (or in general the distribution) 
of such a ratio would be much more useful. To obtain such information from data with 
large dispersion a statistical approach to study the relationship between these two emis- 
sion mechanisms is essential. A statistical study can establish the general trends in the 
data (such as mean values and dispersion) and therefore be helpful in the understanding 
of the general behavior of flares, such as the directionality of the accelerated electrons 
and the magnetic field configuration surrounding the acceleration region. Kane’s (1981) 
investigation dealt with the correlation of the X-ray and type III emissions observed dur- 
ing the solar maximum of the 20 th solar cycle (circa 1969). He found that the correlation 
between X-ray and type III events increased with the intensity and starting frequency of 
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the type III burst, and with the peak energy flux and spectral hardness of the hard X-ray 
burst. He concluded that the observations suggest that the acceleration region is located 
in the corona and that the electrons responsible for both emissions are accelerated in a 
single process. The observed correlations were shown to be consistent with a bivariate 
distribution function of X-ray and type III burst intensities by Petrosian and Leach (1983, 
hereafter PL). Although the distribution could not be determined from Kane’s data, a 
procedure was outlined for the determination of the distribution given a larger data set 
with better resolution of radio burst intensities. Therefore, we decided to look at a large 
sample of flares which had either a hard X-ray event or a type III burst for the most recent 
solar maximum and to apply the method of PL to determine the bivariate distribution 
function and discuss the resulting implications. 

The paper is organized as follows. In the next section, we describe the data set that 
we analyzed. In §111, we use the data to find the resulting bivariate distribution function 
and discuss its properties. A discussion of our results is given in §IV and a brief summary 
in §V. 


II. Data 

The hard X-ray events were obtained from the tables in the Hard X-Ray Burst Spec- 
trometer Event Listing 1980 to 1987 (Dennis et aZ.1988) excluding events which were 
interrupted or were nonsolar in origin. The radio data consisted of observations made 
at Bleien (BLEN) and at the Osservatorio Astronomico Di Trieste (Trieste) tabulated in 
their quarterly report Osservazioni Solari, These data were accumulated and kindly made 
available to us by Dr. A.O. Benz. The BLEN data does not give a precise quantitative 
measure of the strength of the type III bursts such as the flux density. Consequently, we 
use the Trieste data which gives the peak fluxes, the start time, the peak time, and the 
duration of the burst at 237, 327, and 408 MHz for our analysis. 

For our study, a hard X-ray event and a type III burst are considered to be associated 
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if the peak in the radio flux occurred during the duration of the hard X-ray burst or the 
peak in the count rate of the HXRBS detector occurred during the duration of the radio 
burst. We studied bursts which occurred during the maximum of the 21 st solar cycle 
from the beginning of the HXRBS data record 1980 February 19 until 1981 September 
30. During this period 3490 HXRBS events were recorded which are included in our study 
an d we have X-ray spectral information for 241 of these bursts which was kindly made 
available to us by Dr. B.R. Dennis. There were 944 recorded bursts at 237 MHz during 
this same period of time with 134 of these bursts correlated with a hard X-ray burst. 
Table 1 summarizes the statistics of the hard X-ray and the 237 MHz radio data used in 
our analysis below. 


III. Analysis of the Correlation 

a) Distribution of Burst Intensities 

For the purpose of description of the distribution of bursts, we choose the readily 
available peak count rate of the hard X-ray bursts and peak flux density at 237 MHz (in 
solar flux units which is equal to 10 -22 Wm -2 Hz -1 ) for the radio bursts which we denote 
by X and R, respectively. The choice of the 237 MHz channel is somewhat arbitrary, but 
similar results are found for all three radio channels. 

Figure 1 shows the distribution of the 134 correlated events (those with both X-ray 
and radio fluxes) in the R-X plane. This figure shows the wide range of the fluxes in both 
energy bands with no apparent clear correlation between X and R. Instrumental cut-off 
seems to be at 20 or 30 counts per second for X-rays and about 50 to 100 solar flux units 
at 237 MHz. We let the number of bursts with intensities in the range X to X+dX and 
R to R-fdR be Af(X, R)dXdR, where f(X,R) is the distribution function of the bursts 
(normalized to unity) and A is the total rate of occurrence of all bursts. The moments 
of this function are useful for demonstrating the existence of a correlation between the 
occurrence of a hard X-ray burst of intensity X and the occurrence of a type III burst of 
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intensity R. 


From the bivariate distribution f(X , R), the number of hard X-ray bursts with inten- 
sity between X and X+dX with radio intensity greater than R 0 is 

/*00 

N r ,(X) = A f(X, R)dR . (1) 

JRo 

Similarly, the number of type III bursts with intensity between R and R+dR with X-ray 
intensity greater than X 0 is 



f(X,R)dX . 


( 2 ) 


For R 0 = 0, equation (1) gives the distribution N 0 (X) of peak count rates of all X-ray 
bursts. Similarly, N 0 (R ) gives the distribution of radio fluxes of all type III bursts. If the 
radio and the hard X-ray intensities are stochastically independent of each other then the 
distribution is separable; f(X,R ) = Aip(X)(f>(R). For this situation equations (1) and (2) 
become 

roo 

N Ro (X) = AiP(X) <t>(R)dR (3) 

JRo 


and 

N Xo (R) = A<f>(R) r H X)dX , (4) 

JXo 

respectively. In this case the effect of changing the sensitivity R 0 (or X 0 ) is only to 
change the normalization of the distributions, but the functional form is identical to the 
distributions of all bursts N 0 (X ) (or N 0 (R ) ). Therefore, we can immediately test for 
statistical independence by comparing the distribution Nr 0 (X ) or N Xo (R ) for different 
sensitivities R 0 or X 0 with N 0 (X ) or N 0 (R). Figures 2 and 3 show such a comparison 
for the data described in §11. The error bars given are proportional to the square root of 
the number of bursts in each bin (Note that the data is probably incomplete at the low 
end). It is evident that the curves are not parallel and that the magnitude of the slope of 
the distributions N Ro (X) shown in Figure 2 decrease with increasing radio sensitivity R 0 . 
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This means that larger radio bursts are more likely to be accompanied by a hard X-ray 
burst. Alternatively, with the statistically independent assumption, one would expect the 
function G(X/R 0 ) = Nr o (X)/N 0 (X), which is the fraction of hard X-ray bursts which are 
associated with a type III burst of intensity greater than R a , to be constant for all values 
of R 0 . However, as shown in Figure 4, this fraction increases with X indicating that the 
larger hard X-ray bursts are more likely to be accompanied by a type III burst. 

As a qualitative measure of this behavior, we list the best fit power-law indices (ex- 
cluding the incomplete bins) to the curves in Figures 2 and 3 in Table 2. The power-law 
index for N 0 (X ) agrees with the well-known result for the HXRBS data (Dennis 1985 Fig. 
3) and the power-law indices systematically decrease with decreasing sensitivity (increas- 
ing R 0 or X 0 ), both for the hard X-ray and type III distributions. We conclude, therefore, 
that the data is inconsistent with the stochastic independence assumed in equations (3) 
and (4). 

The question that arises is, what kind of correlation exists between these two types of 
emissions or parameters? We have yet to determine the functional form of the distribution 
f(X,R) which gives rise to the observed distributions Nr o (X) and Nx a (R)- In general, 
there axe a multitude of functional forms which can describe a given set of observations. 
Without any loss of generality, we can write 

f(X,R)dXdR = Atp(X)g(Y)dXdY , (5) 

where Y is an arbitrary function of X and R. Here we first analyze the simple form 
Y = X/R discussed by PL. This form for the distribution function gives 

Nr o (X) = A4>(X)G(X/R 0 ) , (6) 

where 

G(Y)= f g(Y')dY' with G(oo) = 1 (7) 

Jo 

is the fraction of bursts with hard X-ray intensity X which have radio bursts of strength 
R > X/Y . Note that N 0 (X) = Aifi(X). This is a useful and general way of describing 
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the correlation between X and R. If X and R are perfectly correlated ( X/R = Y 0 for 
all events) then g(Y ) oc S(Y — Y a ). For complete independence g(Y ) is featureless (e.g. a 
simple power-law). In other words, the “width”of the distribution g(Y) provides a measure 
of the degree of stochastic dependence of the quantities X and R. 

This distribution assumes stochastic independence of the X and Y parameters and 
therefore can be subjected to a similar analysis described above for the X — R distribution. 
The distribution of the correlated events in the Y-X plane are shown in Figure 5. Also 
plotted is the effective instrumental cut-off at R 0 = 75. The absence of data points in the 
upper left region, Y <10 -1 , X > 10 3 of this figure, as well as the absence of data below 
the cut-off line, make these quantities appear to be strongly correlated. However, based on 
the occurrence rates for all flares one only expects one or two points to be in the first of the 
above mentioned regions so that the absence of observed flares there is not unreasonable. 
Due to the presence of a cut-off, the construction of the distributions g{Y) or tj>(X) is 
not straightforward. If X and Y are stochastically independent then the non-parametric 
method described by Petrosian (1986) is one way to construct these distributions. This 
method gives the cumulative functions $(X) = ip(X') dX' and G(Y). In Figure 6 we 
show the cumulative function $(-Y) for the Y parameter in the ranges 0 < Y < 0.4 and 
0.4 < Y < 6.4. The fact that these distributions are fairly similar is an indication of the 
independence of X and Y. The distribution Ny 0 (X) = $(X)G(Y 0 ) versus X is fitted to a 
power-law function from which we calculate the index din ip/d\nX for various values of 
Y 0 . The indices, given in Table 2, show no systematic variation with Y 0 and within the 
errors agree with the power-law “index”found for all bursts ( Y 0 — * oo). This again shows 
that the separation of the distribution given by equation (5) is consistent with the data. 

Using the same non-parametric method we can evaluate the cumulative function G(Y) 
(see Figure 7). The functions ip(X) = N 0 (X)/A and g(Y) = dG(Y)/dY then give a 
complete description of the distribution. This procedure, however, requires an analytic 
form for G{Y). Except for the case of perfect correlation, there is no unique form for 
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g(Y ) or G(F). We chose a few simple functional forms with two free parameters and 
found the values of these parameters that best fit the data. The forms tried were gaussian, 
exponential, and simple power-law like forms. No satisfactory fit was obtained for the 
gaussian or exponential forms. The resulting best fit for the power-law like form is 

G(Y) = 1/(1 + (n/V)‘) . 9 (Y)=bY-\YIY a ) b - l l(\ + (Y/Y 0 ) l f , (8) 

with b = 0.82 ±0.07 and Y 0 = 8 ±2. This fit is shown in Figure 7. To check our method for 
internal consistency, we use the analytic functions g(Y ) and ^(X) oc X~ 1-82 to evaluate 
the expected number of bursts Nr o (X) and compare this distribution with the data. The 
circles in Figure 2 show the results of this comparison for two values of R a . Similarly, we 
can compute the expected distributions Nx„(R ) of type III events, 

N Xo (R) — A f°° ft X)g(X/R)XdX/R 2 . ( 9 ) 

Jx 0 

We see in Figure 3 that the calculated points (circles) agree with the observed histograms 
here as well (The lowest point doesn’t fit due to the incompleteness of the radio data at low 
flux levels.). Therefore, the form for the bivariate distribution function of hard X-ray and 
type III burst intensities given by equation (5) with g{Y ) given by equation (8) describes 
the observed distribution of bursts adequately. 

In order to obtain a more quantitative measure of the degree of correlation between 
the X-ray and type III radiations, we have repeated the same procedure with two other 
simple functional forms for Y: namely, Y' = \/X/R, and Y" = X 2 / R. We find that the 
power-law index of Ny{X ) changes systematically from 1.40 (with 25 points) to 1.63 (with 
132 points) as Y 1 is increased from 0.01 to 1.0. Furthermore, the index of 1.63 does not 
agree with the expected result of 1.82 for all X-ray bursts (cf. Table 2). For the choice 
Y", the power-law index of Ny»{X) is too large with values of 2.44 (28 points), 2.74 (70 
points), and 2.55 (116 points), for Y" = 5, 50, and 10 4 , respectively. Therefore, neither 
Y' nor Y" appear to be stochastically independent of X and our choice of Y = X/R 
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provides a better description of the data. More elaborate functions might be chosen, such 
as Y = ( X — X t h)/R, where X t h may be due to another process contributing to the X- 
ray flux (e.g. thermal bremsstrahlung) but not to the type III emission. In this case, the 
distribution in X t h would also have to be unfolded. We feel that the data does not warrant 
such detailed analysis and that the thermal contribution only adds more dispersion to the 
data (see §IV below). 

We note that our choice for the independent parameters is slightly arbitrary. Instead 
of X and Y we could have chosen R and Y. In fact, if X (or R) and Y were strictly statis- 
tically independent, then one could choose any function of X (or R) and Y as independent 
parameters. This, however, would not alter observable quantities such as the distributions 
Nr o (X), N Xo (R) or the values of the moments of these functions. We have repeated the 
calculations leading to Figures 5 through 7 using R instead of X, and Y as our parameters 
and have fitted the data represented in this way to the function described by equation 
(8) times tp(X) oc X -1-82 = (i?Y) -1-82 and find b = 0.83 and Y 0 = 7. These values are 
in agreement with the values b = 0.82 ± 0.07 and Y 0 = 8 ± 2 found above. It is clear 
from this analysis that the distributions N 0 (X) and N 0 (R ) should be identical. As shown 
in Table 2, the power-law indices (1.82 and 1.75) only differ by slightly more than one 
standard deviation. Therefore, the observed N 0 (R ) distribution could be a random sample 
of the N a (X) distribution. Finally, we note that the ratio of X-ray to radio limiting flux 
Yum = Xii m / Rum ~ 25/75 = 1/3 is less than the representative value Y a of the Y param- 
eter determined by either method. Ideal samples should have Y/, m ~ Y 0 which means that 
we need an improvement of a factor of 24 ( Ru m ~ 3 solar flux units) in the radio data if 
Y 0 ~ 8. Further lowering of Ru m will not increase the size of the data set significantly, or 
improve it, unless accompanied by proportional decreases in Xu m . 

b) Dependence on Spectral Index 

As summarized in Table 1, we have the values of the X-ray spectral index, 7, from 
power-law fits to the X-ray spectra for 241 of the bursts, 34 of which have been observed to 
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be associated with a type III burst. In this subsection, we see to what extent the spectral 
index effects the correlation between the hard X-ray and type III emissions. 


The histogram in Figure 8 shows the fraction of hard X-ray bursts which have detected 
type III fluxes versus the spectral index (solid line) for R > 50. For comparison we also 
show (dashed line) the data from Figure 10 of Kane (1981). Both histograms show that 
a larger fraction of harder (smaller 7) X-ray bursts produce type III bursts. This effect is 
much less pronounced in the new data. In fact, a x 2 -test gives x 2 = 4.6 (with 7 degrees of 
freedom) and therefore a probability of 0.70 that the x 2 found from a parent distribution 
which is uniform would be this large or larger. Hence, we cannot exclude the possibility 
that there is not a correlation. This presumably is not true for Kane’s (1981) results and 
we do not know the reason for this difference. To determine whether or not the spectral 
index dependence that we do see is due to selection effects, such as correlation of X-ray 
intensity with spectral index, we have preformed the following test. We binned all X-ray 
bursts with known spectral index (241 of them) by spectral index. Since the probability 
that an X-ray event of strength X is associated with a type III burst is known (cf. Figure 
2), we can compute the expected number of type III associated events for each bin. This 
is found by summing the probability of type III association for each individual event over 
all X-ray events in the bin. The result of this calculation is shown in Figure 9 along 
with the observed number of hard X-ray bursts (with known spectral indices) associated 
with a type III burst. We see that for 7 > 4.5 the number of associated bursts is equal 
to that expected from the total distribution while for 7 < 4.5 the number of associated 
bursts is twice that expected. In this case, we find x 2 = 3.0 (8 degrees of freedom) which 
gives a higher probability (0.93) that this x 2 could have come from a distribution which 
is uniform in spectral index. This indicates that some sampling bias exists. In order to 
conclude that the increased correlation between hard X-rays and type III emission for 
harder X-ray spectra is a real effect, a larger sample is needed. 
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IV. Discussion 


We have shown that the observed distribution of the type III and hard X-ray emission 
from solar flares can be described by a bivariate distribution with independent parameters 
the X-ray intensity and the ratio of the X-ray to type III intensities. This implies a 
correlation between these two types of emission. However, this correlation is stochastic 
as demonstrated by the large range (or “width”) of the distribution g(Y) in equation (8). 
Strictly speaking, both the average value and mean dispersion of Y for this distribution 
diverges. This, however, is an artifact of the simple distribution chosen here and can be 
remedied by the introduction of additional parameters. The important conclusion we draw 
from this is that a typical value for the ratio Y is the quantity Y 0 ~ 8 and the range of Y 
is even greater than Y a and can be as large as 10 2 F o ~ 10 3 . We now discuss the meaning 
of these results and the results from consideration of the spectral index dependence in the 
frame work of models for the impulsive phase of the flare and for the two types of radiation 
under consideration here. 


a) Nature of the Correlation 

We first consider the correlation we have found between the X-ray and type III emis- 
sions and in particular the appropriateness of the choice Y = X/R as a stochastic param- 
eter. 

The X-rays of energies 20 ~ 100 keV are produced mainly by bremsstrahlung of 
electrons with comparable energies. Three distinct models have been advanced to explain 
the solar hard X-ray emission by electrons, the thick-target, thin-target, and thermal 
models (trap plus precipitation models are classified as thick-target here because the X- 
ray production is via thick-target bremsstrahlung of the precipitating electrons). In both 
the thick and thin target models, nonthermal electrons produce the impulsive X-ray bursts 
and the strength of the burst is proportional to, n x , the number of nonthermal electrons. 
In thermal models, the level of X-rays produced is proportional to n 2 . Since the corona is 
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optically thin to photons in this energy range, the observed radiation is also proportional to 
n x (nonthermal) or n 2 x (thermal). Although there is often some contribution from thermal 
emission, the impulsive X-ray bursts are most likely due to nonthermal electrons. 

The relationship between the number of type III producing electrons, n///, and the 
intensity of type III bursts is more complicated. Type III bursts are produced by nonther- 
mal electrons with energies E ~ 30 keV ()3 ~ 1/3) which develop a positive slope in their 
velocity distribution as they propagate. This situation is unstable to the generation of lon- 
gitudinal plasma waves (Langmuir waves) with frequency near the local plasma frequency. 
These waves don’t propagate or escape the plasma, so that a mechanism for conversion of 
these waves to escaping radiation is needed. There are two processes which are believed to 
occur. The coalescence of two longitudinal waves to form a transverse wave at twice the 
plasma frequency (second harmonic emission) and the scattering of a longitudinal wave 
off ion-acoustic turbulence into a transverse wave with frequency near the local plasma 
frequency (fundamental emission). To find the relation between n/// and the intensity of 
type III radiation we must first relate n/n to the level of plasma turbulence generated, 
Wi, and then relate this to the intensity, R, of type III radiation. 

The first of these relations is simple Wi oc n/ 77 . The relationship between the tur- 
bulence level and the intensity of type III radiation was discussed by Smith and Spicer 
(1979). For second harmonic emission they showed that R oc Wi when optically thick and 
R oc Wf for optically thin sources. Therefore, R oc n/// or R oc njjj for second harmonic 
emission from an optically thick or optically thin source, respectively. They also showed 
that R oc Wi for optically thin fundamental emission. For optically thick sources the 
fundamental emission depends on the ratio of Langmuir turbulence to ion-acoustic turbu- 
lence. However, the second harmonic emission is usually dominant for type III bursts and 
for solar type III bursts the source is generally optically thick. So that for most type III 
bursts R oc W[ oc n///. 

Combining the results for X-ray and type III emission, we see that for the most 


13 



likely situation where impulsive hard X-rays are produced by nonthermal electrons and 
the type III burst is second harmonic emission from an optically thick source the quantity 
Y = X/R oc n x /njji. Therefore, if this ratio is independent of X or R then Y is a good 
variable with which to characterize the correlation of X-ray and type III bursts. Other 
situations are of course possible and would lead to other parameterizations. For example, 
for nonthermal models (X a n x ) with optically thin second harmonic type III emission 
(. R oc n 2 in ) parameterization in terms of Y" = X 2 /R oc ( n x /nm ) 2 would be more natural. 
Or, if X-rays are produced by thermal bremsstrahlung {X oc. n 2 ) and type III radiation 
by optically thick second harmonic emission, then the parameter Y 1 = y/X / R oc n x /nuj 
would be the natural parameter if nm oc n x . For these cases, however, there may not be 
any correlation between n x and nm so that one would expect no correlation between X 
and R and a distribution function as described in connection with equations (3) and (4). 
As shown in the previous section, all of these other possibilities are less likely than the 
parameterization in terms of X and Y = X/R, corresponding to nonthermal X-rays and 
optically thick second harmonic type III emission. 

Thus, we can conclude that the ratio Y = X/R is proportional to the ratio Y e = 
nx/nm of the number of X-ray to type III producing electrons; Y = rjY e . Part of the 
large (three or four orders of magnitude) observed dispersion in the value of Y is due to 
the dispersion of the proportionality “constant ” 77 which depends on the properties of the 
flare plasma, especially those determining the level of plasma turbulence W\. We estimate 
rj to be 10 -1 to 10 -3 based on the relation between n x and X which is straightforward 
(cf. Hoyng, Brown, and Van Beelc 1976) and on the relation between nm and R which 
is very model dependent. We have used the relations given by Smith and Spicer (1979) 
and assumed a velocity of c/3 for the electrons, a range for plasma temperature of 10 6 
to 10 7 , and efficiency of plasma wave generation of 0.1 to 1 . Equally important in the 
dispersion of Y is the variation of the ratio Y e from flare to flare. The ratio Y e depends 
on a number of factors related to the acceleration mechanism and the transport of the 
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accelerated electrons in the solar atmosphere. The most important among these factors 
is the ratio of the number of accelerated electrons which end up on open field lines (lines 
extending to the upper corona and/or the interplanetary medium) to those on closed field 
lines, which we denote by e. 

All electrons on the closed field lines and a fraction f x of those on the open field lines 
contribute to the X-rays while the rest contribute to the type III emission. (We assume 
that X-ray emitting electrons have a negligible contribution to the type III emission and 
vice versa). Thus, the ratio Y e = (1 + ef x )/e( 1 — f x ). The ratio f x in turn depends 
on two factors. The first of which is the fraction of the electrons on the open field lines 
directed outward from the sun, f up , which is determined by the pitch angle distribution 
of the accelerated electrons (e.g. f up = 1/2 for isotropic pitch angle distribution). Due to 
pitch angle diffusion and field convergence a fraction f re f of the down going electrons are 
reflected back so that f x — (1 — / re /)(l — f up ) and 


1 + e(l - f ref )( 1 - f up ) 

^(fup d" fref fref fup) 


( 10 ) 


which is the same as equation (17) of PL with different notation. 

First we note that based on our estimated range of values for 77 , a typical observed 
ratio Y = Y 0 requires the value Y eo ~ Y 0 jr\ ~ 10 3±1 . The important point here is that 
despite the large uncertainty we require Y e 1 even for those flares with minimal X-ray 
emission having Y ~ 1 . According to equation (10), this can be the case if either e < 1 or 
(fup + fref) < 1 or both. It can be shown that within a factor of two this expression can 
be simplified as 


I g 


1 + 6 


t(fref d” fup ) 


( 11 ) 


which should be sufficient for our discussion considering the large uncertainties in the other 
parameters. The large dispersion in the value of Y can then be easily accommodated by 
the dispersion in the three parameters in equation (11). Of these, e and f up depend on 
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the acceleration process and the geometry of the magnetic field (0 < e < oo, 0 < f up < 1) 
while f re f depends on the transport effects and is expected to be small (f re f < 1). 

Another contribution to outgoing type III producing electrons may come from possible 
diffusion across the magnetic field lines (from closed to open field lines) as in the model 
proposed by Emslie and Vlahos (1980, hereafter EV). In the general version of the model, 
Y e is given by equation (10) or (11) with the addition in the denominator of another term 
fdiff which is the fraction of electrons in the closed loop that diffuse to the open field lines 
with pitch angles so that they stream outward from the sun. In the version of this model 
proposed by EV, all particles are on closed field lines so that e = 0 and 

n * ■ (i2) 

In general, fdiff is expected to be small so that Y e 1 as required, but the distribution 
of fdiff (the values of the expected average and rms values) is not explicitly specified by 
the model and difficult to predict. 

Finally, let us consider a third model for the production of type III radiation pro- 
posed by Sprangle and Vlahos (1983, hereafter SV). In this model, the primary electrons 
which are accelerated in a closed loop with a converging magnetic field reflect to acquire 
a loss cone distribution. This electron distribution is unstable and results in cyclotron 
maser emission. The emitted radiation then escapes the closed field lines of the loop and 
accelerates secondary electrons on open field lines. This model is able to account for short 
time delays between peaks in the type III emission relative to peaks in the X-ray emission 
(Dennis et al. 1984). SV show that the number of accelerated electrons is proportional to 
the magnitude of the loss cone, and therefore to the number of reflected electrons, times 
factors related to the efficiency of the secondary acceleration process. Since the number of 
reflected electrons is proportional to n x , this model also agrees with the Y parameteriza- 
tion that we used and predicts a value of Y e given by equation (10) or (11) with addition 
of another term a f re f in the denominator for the contribution to n/// by the secondary 
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process. Of course, as in the EV model, this model is of interest when e = 0. So that in 
the simplest form of the model Y e oc As we shall see below f re f <C 1 and Y e 1 as 

expected. However, again determination of the average or rms value of Y e is difficult in 
this model. 

b) Dependence on Spectral Index 

The discussion above shows that most model predictions agree with the observed 
distribution and its parameterization in terms of X and Y = X/R. We now consider the 
dependence of the X-ray and type III correlation on the X-ray spectral index. In general, 
for nonthermal models there is a simple relation between the X-ray spectral index 7 and 
the electron spectral index S. In particular, for the thick target models considered here in 
their simplest form, 7 = £ — 1. This relation is not exact and depends on the parameters 
of the model (see Leach and Petrosian 1983, or Leach 1984), but for our purposes here it is 
sufficient to note that as 8 increases 7 increases so that the observed correlation between 
the distribution and 7 (if any) reflects correlation between 8 and Y e . 

As shown in the previous section, the new data does not show the strong dependence 
on spectral index found by Kane (1981), so that the degree of the correlation is not well 
determined. However, we believe that bursts with harder X-ray spectra are more likely to 
also have type III emission. This can be described mathematically by making Y 0 (or b ) to 
be a function of spectral index. If we take b to be independent of spectral index then for 
each spectral index bin the fraction of bursts with type III emission is 

fgijX)Gl X/RM^dX 
/“ i,(X) iX 

For ip(X) oc X~ a and G(Y) given by equation (8) this reduces to 

r> OO 

f(i) = zr'Ro, - 1) / dZ Z‘- 7(1 + Z‘) , (14) 

Jz 0 

where Z 0 = Yu m /Y 0 ( 7). As noted above Yu m ~ 1/3 <C Y 0 so that Z 0 <C 1, in which 
case, to within some slowly varying logarithmic factor, the variation with spectral index 
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of this fraction can be approximated as /( 7) \Yu rn /Y 0 ('y)] b . Our results and those by 
Kane (1981) shown in Figure 8 indicate /(7) « 7- 1 and 7 -3 , respectively. If we assume 
/( 7) « 7 -2 and neglect the slight difference between the exponent b and unity the data 
would require F 0 (7) oc 7 2 . This is a very approximate relation that implies Y e oc 6 2 , a 
requirement that the models must satisfy. If instead of Y 0 we assume b to be a function 
of 7 then the above approximation is still valid and In f(j) = b(~f)ln(Yu m /Y 0 ). Thus, an 
observed /("/) cx y~ a can be described by b = b 0 — aln(7/7 0 ), where b 0 is the value of b at 
7 0 . In what follows we consider the case that b is constant and Y 0 = Y 0 { 7). 

In the simple nonthermal model it is evident from equation (11) that we require e ^ 0 
otherwise there will not be any type III emission. Furthermore, in its simplest form this 
model does not require or predict dependence on the spectral index 8 for this parameter 
or for f U p which is a moment of the pitch angle distribution of the accelerated electrons. 
Then, as proposed by PL, the only dependence on 8 will arise from the variation of f re f 
with 6. As also shown by PL, there is considerable variation of / re / with 8 especially 
for models with electrons beamed primarily downward f up — * 0 so that Y e oc From 

Table I of PL we see that for this case (their parameter oil < 1) fref increases more rapidly 
than linearly, so that the trend agrees roughly with the above requirement. 

On the other hand, if the pitch angle distribution of the accelerated electrons is 
isotropic f up = 1/2, then we require e < 1 so that Y e oc (1 + 2/ re /) _1 . Since, in gen- 
eral f re f < 1, this predicts a slower variation. In addition, for an isotropic pitch angle 
distribution the above mentioned Table of PL shows a slower variation of f re f with 8. 
Consequently, we expect very little variation of Y e with 8 and Y 0 (or /(q) ) with 7. Note 
that for both of these cases the dependence on 8 is weakened with increasing magnetic 
field convergence. 

For a quantitative comparison of these models with the data we assume that the 
spectral indices are related as in the simple thick-target model, 7 = 6 — 1, and choose e/77 
so that Y 0 = 8 for the median spectral index 7 = 4. This means that r/Y e (8 = 5) = 8. 
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The fraction /( 7 ) is computed using equation (13) with Y 0 ( 7 ) given by equation (10). The 
model parameters are listed in Table 3 and the corresponding /( 7 )’s are shown in Figure 
10 for comparison with the data. We find that the calculated fractions agree with the 
above qualitative discussion. The isotropic models show very little variation with spectral 
index while the beam models fit the data well, suggesting that the pitch angle distribution 
of accelerated electrons is beamed toward the sun. 

In the EV model (with e = 0) Y e oc /Ji//- It is not clear why fdiff would be dependent 
on spectral index 8. However, since higher energy particles penetrate to higher depths 
below the transition region and if the diffusion rate across the field lines also increases 
with depth then bursts with harder spectra (lower 8) will have larger fdiff , a lower Y e or 
Y 0 , and therefore a larger fraction /( 7 ) associated with type III emission. 

In the SV model, the strength of the loss cone determines the strength of the type III 
emission. Since the bulk of the strength of the loss cone comes from lower energy electrons 
and since the loss cone is produced by mirroring before much collisional energy loss occurs, 
we expect little dependence on the spectral index. For this model the field convergence is 
strong so that f re f approaches unity with very weak dependence on 8. 

Clearly we require more spectral index information before such observations can be 
used to distinguish between the various possibilities described above. 

V. Summary 

We have analyzed a large number of hard X-ray events which occurred during the 
maximum of the 21 s * solar cycle in order to study their correlation with type III bursts. 
From our analysis, we found that the distribution of occurrences of hard X-ray bursts which 
are correlated with a type III radio burst is significantly different than the distribution of 
all hard X-ray bursts, which is consistent with results found by previous authors. We 
were able to show that this result is consistent with the assumption that the hard X-ray 
and type III intensities are somewhat correlated. A bivariate distribution function of the 
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burst intensities was fit to the data. From this distribution we determined that the typical 
ratio of X-ray intensity to type III intensity is Y ~ 8 and that the ratio of the number 
of X-ray producing electrons to type III producing electrons is in the range Y e ~ 10 2 to 
10 4 . However, there is a large dispersion in the ratio of peak X-ray intensity to type III 
intensity. 

Three models that have been proposed to explain the relation between hard X-ray 
and type III producing electrons were examined in the context of these observations. All 
models predict that Y oc Y e and that it is independent of both X-ray and type III intensity 
which agrees with the data. The important differences in the models are whether or 
not both electron populations are accelerated simultaneously and whether the correlation 
between hard X-rays and type III radiation increases for harder X-ray bursts. The data 
indicate that this correlation increases slightly with decreasing X-ray spectral index. The 
increased correlation of X-ray and type III events with decreasing photon spectral index 
is explained in the model of PL by the increase in the fraction of reflected electrons for 
harder electron spectra suggesting that reflected electrons play a role in type III bursts. 
The observed variation with spectral index is best fit by models in which the accelerated 
electrons, especially those on open field lines, are beamed toward the sun. More data is 
needed to further quantify this correlation so that the predictions of the impulsive phase 
models can be further tested. However, this also requires more detailed acceleration models 
so that the spectral index dependence of the ratio Y t can be determined. 
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Table 1 


Summary of hard X-ray, type III data, and the correlated events. 

Hard X-Ray Bursts 


All events: 

Total number of bursts 3490 

Number of bursts with known spectral index 241 

Correlated events: 

Total number of bursts 134 

Number of bursts with known spectral index 34 

Type III Bursts 

Total number of bursts 944 

Number of correlated bursts 134 
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Table 2 


Power-law fit to the functions Nr o (X), Nx a (R )> and Ny o (X). 


Nr.(X) N x .(R) Ny.(X) 


R 0 

n 

a 

A a 

X 0 

n 

a 

Aa 

Y 0 

n 

a 

Act 

0 

3490 

1.82 

0.02 

0 

944 

1.75 

0.04 

oo 

3490 

1.82 

0.02 

75 

126 

1.39 

0.06 

20 

132 

1.54 

0.08 

6.4 

113 

1.82 

0.07 

150 

93 

1.33 

0.06 

40 

101 

1.48 

0.09 

1.6 

o 
o 
1— 1 

1.78 

0.08 

300 

63 

1.30 

0.07 

80 

76 

1.44 

0.10 

0.4 

59 

1.98 

0.17 

600 

38 

1.27 

0.09 

320 

36 

1.37 

0.13 

0.1 

24 

1.95 

0.25 


n is the number of points in the distribution. 
a is the power-law index such that f(z ) cx z~ 
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Table 3 


PL Model Parameters 


model 

fup 

r m 

log(e/77) 

51 

0 

1 

2.90 

B5 

0 

5 

0.37 

11 

1 

2 

1 

-0.62 

15 

1 

2 

5 

-0.67 


r m is the factor by which the magnetic field increases from the injection point to the 
transition region. 
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Figure Captions 


Fig. 1. Distribution of the correlated events in the R-X plane. 

Fig. 2. Distribution of all HXRBS events N a (X ) and Nr 0 (X) for R a = 75 and 300 with the 
R 0 = 300 histogram reduced by a factor of ten for clarity. The solid lines show the 
best fit power-law for each of the histograms. The circles are the values calculated 
using the best fit to the function G(Y). 

Fig. 3. Same as Figure 2 except for type III distributions N 0 (R ) and Nx 0 (R) for X Q = 40 
and 320 with the X Q = 320 histogram reduced by a factor of ten for clarity. 

Fig. 4. Fraction of X-ray events correlated with a type III event for minimum radio 
intensities of R a = 75 and 300 solar flux units. The R„ = 300 histogram is reduced 
by a factor of ten for clarity. 

Fig. 5. Distribution of the correlated events in the Y-X plane. Also shown is the instru- 
mental cut-off at R=75. 

Fig. 6 . Cumulative distribution $(AT) for 0 < Y < 0.4 and 0.4 < Y < 6.4 obtained using 
the non-parametric method with the 0 < Y <0.4 curve reduced by a factor of two 
for clarity. Also shown is the best fit power-law for each of the distributions. 

Fig. 7. The cumulative distribution, G(Y), obtained using the non-parametric method 
along with the best fit curve equation ( 8 ). 

Fig. 8 . Ratio of correlated events to total HXRBS events as a function of photon spectral 
index (solid). Also depicted are the results from Figure 10 of Kane (1981) (dashed). 
The power-law functions / oc 7 -1 and are also shown. 

Fig. 9. The number of observed correlated bursts (solid) and the expected number of 
correlated bursts (dashed) as a function of X-ray spectral index. 

Fig. 10. The observed fraction of correlated events as a function of spectral index /( 7 ) is 
shown (histogram) along with the predictions of the PL models of Table 3. Models 
B 1 and II are shown as solid lines while B5 and 15 are dashed. 
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Count i n g Rd t e 


X) and Nji o (X ) for R„ — 75 and 300 with the 
>r of ten for clarity. The solid lines show the 
>grams. The circles are the values calculated 




Fig. 3. Same as Figure 2 except for type III distributions N„(R ) and Nx„(R) for X a = 40 
and 320 with the X a = 320 histogram reduced by a factor of ten for clarity. 




Peak X-Ray Counting Rate 


Fig. 4. Fraction of X-ray events correlated with a type III event for minimum radio 
intensities of R 0 = 75 and 300 solar flux units. The R 0 = 300 histogram is reduced 
by a factor of ten for clarity. 









.5 3.0 3.5 4.0 4.5 5.0 5.5 6.0 6.5 

Photon Spectral Index 

Fig. 8 . Ratio of correlated events to total HXRDS events as a function of photon spectral 
index (solid). Also depicted are the results from Figure 10 of Kane (1981) (dashed). 
The power-law functions / oc 7 “* and 7~ 3 are also shown. 



Number of Bursts 




Fraction Correlated 



Fig. 10. The observed fraction of correlated events as a function of spectral index /( 7 ) is 
shown (histogram) along with the predictions of the PL models of Table 3. Models 
B 1 and II are shown as solid lines while B5 and 15 are dashed. 





